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Resume An exhaustive kinetic model for the atoms in a ID Magneto-Optical Trap is derived, without 
any approximations. It is shown that the atomic density is described by a Vlasov-Fokker-Planck equation, 
coupled with two simple differential equations describing the trap beam propagation. The analogy of such 
a system with plasmas is discussed. This set of equations is then simplified through some approximations, 
and it is shown that corrective terms have to be added to the models usually used in this context. 
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1 Introduction 

The magneto-optical trap (MOT) is the primary tool 
to cool atoms. The development of this technique led to 
spectacular breakthroughs in experimental quantum phy- 
sics : MOTs are the first step in the realization of optical 
lattices [T], cold molecules [5], or Rose- Einstein conden- 
sates [3]- But the MOT is also an interesting object per 
se. It produces a cloud of cold atoms, the physics of which 
is complex. In particular, spatio-temporal instabilities of 
this cloud are commonly observed |4I5) . Several models 
with different approaches have been proposed to describe 
these dynamics, and to identify the mechanisms leading 
to instabilities. Unfortunately, none of these models gives 
a satisfying description of the observed dynamics. In [|], 
a very simple model allowed to describe experimentally 
observed instabilities. This model, improved in [B], pre- 
dicted another type of instabilities, which were effectively 
observed in |7l8j . However, the agreement between the mo- 
del and the experiments was only qualitative. Moreover, 
it concerned the particular case where the counterpropa- 
gating trap beams are obtained by retro-reflection, i.e. a 
global asymmetry is introduced in the trap. 

Recently, it has been proposed to describe the sym- 
metric MOT, with all trap beams which are independent, 
as a weakly damped plasma 5 9 . Indeed, the cloud of 
cold atoms in a MOT is a confined dilute object with 
long-range interactions, as in plasmas |10lll| . This model 
predicted the existence of instabilities above a relatively 
high threshold, so that instabilities should exist only in 
large MOTs. This seems in contradiction with the obser- 
vations related in p]. Moreover, no direct comparison with 
the experimental temporal regimes allowed validating this 
model. More recently, a more complete description was 
derived using the methods of waves and oscillations in 
plasmas, leading to interesting predictions [T^. But as in 
the previous case, this study is based upon intermediate 



well-established results, valid only in specific cases (e.g. a 
low beam intensity or a negligible viscosity) . These condi- 
tions do not correspond in general to the experimental 
situations, and indeed, these results were not compared 
to experimental results. 

It appears from these numerous works that a reference 
model for MOT atom clouds lacks. Such a model should 
be as general as possible, and should at least describe the 
usual experimental situations. In particular, such approxi- 
mations as the low saturation limit should not be done a 
'priori. Another interest of such a model is that it could 
help in determining precisely the analogies between MOT 
atom clouds and other systems, such as plasmas. 

The present work is a first step towards such a model. 
Its aim is to build a model, with the least possible hypo- 
theses and approximations, of a ID symmetric MOT. The 
resulting set of equations describes as precisely as possible 
the dynamics of atoms inside a ID MOT, and constitutes 
a basis model. If simplifications are necessary, approxima- 
tions should be applied on this model, a 'posteriori. 

The paper is organized as follows : section 2 defines 
the bases of the model, while in section 3, we show that 
the dynamics of the MOT phase space density can be des- 
cribed by a Vlasov-Fokker-Planck equation with two re- 
laxation processes. In section 4, we derive the different 
terms of this equation of evolution as a function of the 
usual experimental parameters. In section 5, we establish 
the equations of propagation of the trapping beams. We 
obtain a set of coupled equations fully describing our sys- 
tem, in the general case. This set of equation can be solved 
numerically to obtain solutions of the general system. Ho- 
wever, to have a better understanding of the MOT, it is 
interesting to go further in the analytical approach, sim- 
plifying a posteriori the model, as discussed above. This 
is done in section 6, where we examine an approximation 
for the atomic response. 
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Figure 1. Sketch of the ID system considered here, a) the 
two contrapropagating beams have opposite circular po- 
larizations, b) The laser beams interact with "three-level" 
atoms onaJ = 0— >-J' = l transition, the degeneracy of 
which is lifted by an inhomogeneous magnetic field. Q± 
are the Rabi frequencies associated with the beams ; A± 
are the effective detunings (see section 4.1). 



2 Definition of the model 

As discussed above, we consider here a ID configura- 
tion. Two counterpropagating laser beams with opposite 
circular polarizations interact with the atoms, as sketched 
on Figure [1^. The beam with the a^ polarization comes 
from the negative abscissa, and is denoted by the minus 
sign (intensity I-{x, i), wave vector k-). In the same way, 
the beam with the a~^ polarization comes from the positive 
abscissa (intensity I+{x,t), wave vector fc+). Forces origi- 
nate from the exchange of momentum between the atoms 
and the electromagnetic field. We consider here that the 
atoms are the simplest ones for which the magneto-optical 
trapping is possible. The laser frequency wl is tuned in the 
vicinity of a ./ = -^ J' = 1 transition with a frequency 
Wo (Fig. [lb). 

In the 2D phase space, a point has the coordinates 
{x,p), where x is the position and p is the momentum. To 
describe the cloud of cold atoms, we introduce the phase 
space density p (x, p, t). Formally, the complete atomic sys- 
tem is described by a phase space density including both 
the hot and the cold atoms. However, laser cooling acts 
only in a limited region of space, typically the intersec- 
tion of the laser beams, and for moderate velocities. As a 
consequence, the surrounding hot atoms are considered as 
a large reservoir which remain in a thermal equilibrium at 
room temperature. In the following, we neglect the action 
of the cold atoms on the hot background. These approxi- 
mations lead to a huge simplification in the description 
of the collisions. Thus, to describe the cold atom dyna- 
mics, we have just to derive the equation of evolution of 



p {x,p, t) from the basic principles of atomic physics. This 
approach allows us to go beyond the Wieman description 
of multiple scattering in terms of absorption cross sec- 
tions |10I13] . p{x,p,t) influences the propagation of the 
trapping beams, and reciprocally the modification in the 
beam intensities changes the evolution of the phase space 
density. Thus, we expect to obtain a system of coupled 
nonlinear differential equations. As our aim is to build a 
theoretical frame in which most of the experimental si- 
tuations can be explored, we limit as much as possible the 
initial approximations. In particular, we do not restrict 
our study to the low saturation limit, as in ]TI\ . 

To derive the system of coupled equations, we proceed 
in two steps. We first assume that we know the intensities 
of the beams everywhere in the sample, and we derive the 
equation of evolution of the density in phase space (section 
3 et 4). Then, we write the equations of propagation of the 
beams assuming that we know the atomic density in phase 
space (section 5). 



3 Evolution of the phase space density 

To derive the equation of evolution of the phase space 
density p {x,p, t), we consider an elementary cell centered 
in (xo,po), with dimensions Sx and Sp. The number of 
atoms contained in this cell is N{t) = p{xo,po,t) SxSp, 
where we assume that Sx and Sp have been chosen small 
enough to neglect higher order corrections. The variations 
of N between t and t + St are governed by three distinct 
phenomena : (i) an atom is kicked in or out by a collision 
with the hot background gas or (ii) an atom crosses the 
border of the cell because of the evolution of either its 
position or (iii) its momentum. The simultaneous change 
in position and in momentum would lead to second order 
terms in St, which are neglected. 

The coUisional processes lead to two terms, one for 
losses and one for gains : 



NSt , , , , 
SNcoii = h A da; op or 



(1) 



where r is the mean time interval between two collisions 
with the hot atoms. The source term A is due to the col- 
lisions between two hot atoms that tend to restore the 
Maxwell-Boltzmann distribution for the phase space den- 
sity (uniform in x, gaussian in p). A contribution from 
the collisions between cold atoms could also be considered 
but, in the conditions where instabilities are observed in 
a MOT, this last contribution can be neglected. 

The second mechanism is the variation of N due to the 
velocity of the atoms. As depicted in Fig. [2^, the atoms 
that cross the borders in xq ± Sx/2 with a velocity po/m, 
where m is the atomic mass, between t and t + St are those 
in the close vicinity of that border. The area of these zones 
is simply SppoSt/m. The variation SN^ in N due to the 
velocity of the atoms can then be written : 

.,. SppoSt r / Sx \ / Sx , 

t>iVx = p I a;o - — ,po,M -p I a;o + Y'^O'^ 

= -{SxSpSt) —--p{xo,po,t) (2) 

m ox 
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Figure 2. Modification of the number of atoms contained 
in the elementary cell due to the evolution of (a) the po- 
sition and (b) the momentum. The dashed regions contri- 
bute to losses, while dotted regions contribute to gains. 
In (a), all the atoms in the indicated regions will cross 
a boundary (we assume that po > 0), and the areas are 
proportional to St. In (b), a small fraction of the atoms, 
proportional to 6t, will receive a momentum kick. We have 
to consider separately the momentum kicks with q > and 
those with q < 0. 



The last mechanism that changes N is due to the 
changes in atomic momentum, which is more tricky to 
evaluate. The atoms undergo cycles where one photon is 
absorbed and another one is emitted. We consider that all 
the underlying physics is described by a probability per 
unit of time V {x,p^t,q) : the probability to change the 
atomic momentum by a quantity q between t and t + St is 
simply StV {x,p, t, q). We assume that the time interval St 
is small enough so that an atom can undergo at most one 
photon scattering event. In such a case, q = h{ka — kg) 
where ka and kg are respectively the wavevector of the ab- 
sorbed and emitted photon. As a consequence, q is boun- 
ded to the range [—2hk,2hk]. A fraction, proportional to 
St, of the atoms contained in this region will change their 
momentum by q. On the side p+ = po + ^p/2, the atom 
number variation SNp^ is the difference between the in- 
coming atoms {q < 0) and the lost atoms {q > 0) : 



SNp+ ^ [SxSt) / dq dq'{Af+ -A/'_) (3) 
Jo Jq 

with A/± — V {x,p+ ± q',t, ^q) p{x,p+ ± q',t) 



To simplify this expression, we assume that the product 
V {x,p, t, q) p {x,p, t) varies slowly with p. Then, the inte- 
grand can be approximated by its Taylor expansion in the 
vicinity of p+. To recover the momentum diffusion pro- 
cess responsible for the non-zero temperature of the trap- 
ped atoms, we have to expand the integrand up to first 
order in q' . The inner integration is then straightforward 
because the integrand becomes a linear function of q' . The 
same calculation has to be done on the opposite side in 
p_ = Pq — Sp/2, which gives an analogous expression for 
SN^.. 



It appears natural to introduce the mean values of q 
and q^ : 



/2hk 
qdqV{x,p,t,q) 
-2hk 



(4) 



-2hk 
2hk 



D{x,p,t)= i - dqV{x,p,t,q) (5) 

J-2hk ^ 

These quantities can be interpreted as the mean force and 
the momentum diffusion coefficient. The variation SNp — 
(5iVp_|_ -|- SNp- in N due to the change in momentum can 
be written as : 



SNp = -{SxSpSt) 



- (FP) - 1^ (Dp) 



dp 



dp^ 



(6) 



Collecting all the contributions ([T]) , ^ and ([5]) and going 
to the limit St ^ 0, we obtain the equation of evolution 
for the atomic density 



d p d d d 

-.P+-—P+ — iFp)-^{Dp) 



dt 



m dx 



dp 



dp^ 



T 



(7) 



The three left terms of this equation are characteristic 
of a Vlasov type kinetic model, except for the third term, 
as here the force depends on the velocity of the particles. 
Thus this last term is a drift term, and together with the 
fourth term, which is a relaxation term, they denote a 
Fokker-Planck description. So the motion of cold atoms 
in a MOT appears to be described by a Vlasov-Fokker- 
Planck (VFP) equation. Thus MOTs are part of a large 
class of systems described by the VFP equations, as e.g. 
plasmas [14] , stars [15] , granular media [16] or electrons in 
a storage ring [T7] . The different systems are characterized 
by the dependence of the force F on the phase space den- 
sity p. For example, for a plasma without magnetic fields, 
the Vlasov-Poisson- Fokker-Planck is used. 

In these systems, the Fokker-Planck terms denote the 
immersion of the particles in a thermal bath, which corres- 
ponds for the cold atoms to the laser light. But in the case 
of cold atoms, a second bath, namely the hot atoms of the 
residual gas, produces a second relaxation term, as well 
as a source term (the right hand terms of Eq. [7]) : indeed, 
the standard MOT is an open system, where the total po- 
pulation can vary. This last point is crucial and cannot be 
neglected, as the collision processes between hot atoms are 
the one and only source of velocity redistribution allowing 
for a high MOT population. 

Because of the two relaxation processes involved in the 
MOTs, the question of the mechanism in which originate 
the instabilities is far from being trivial. Indeed, it is well 
known that the VFP equation may have, for adequate 
parameters, unstable solutions. But as in JTj, instabilities 
lead to variations of population, it is clear that the loss and 
source terms can also generate unstable solutions. There- 
fore, instabilities may originate in two different mecha- 
nisms, and it would be interesting to search if one of them 
prevails. An interesting case is that of the lower MOT in 
a double cell system. In this case, the gas pressure in the 
cell is so low that the collision losses are negligible. If the 
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loading process by the upper MOT is stopped, Eq. (O 
becomes fully equivalent to plasmas equations, and only 
one relaxation process remains. Thus it would be inter- 
esting to study experimentally such a system, as the fact 
that instabilities persist or not in this case could help in 
understanding the origin of these instabilities. 

To go further, we need to evaluate the force Q and the 
diffusion coefficient ([5]) for a specific atom. To do so, we 
have to calculate the probability V {x,p,t,q) associated 
with the momentum exchange q = h{ka — fcs)- However, 
as the probability to emit a spontaneous photon in the 
direction fc^ is always identical to that for the opposite 
direction —kg, the contribution of spontaneous emission 
cancels in F, while in D the cross-term in q^ vanishes. As 
a consequence, it is not necessary to evaluate completely 
V {x,p, i, q), but only the absorption rate for a photon of 
wavevector ka- 



4 Radiative forces and coefFicients of diffusion 

As stated above, we consider the simple scheme of a la- 
ser frequency wl tuned in the vicinity ofaJ = 0— J-J' — 1 
transition with a frequency ujq. These atoms are not sui- 
ted for sub-Doppler cooling mechanisms [18], but these 
subtle processes are not needed to describe the tempera- 
tures measured experimentally in the cold atomic cloud 
obtained with high intensities and small detunings for the 
trapping beams [7j . In the case of retroreflected trapping 
beams, this simple atom allows to develop a model which 
convincingly reproduces the experimental observations of 
instabilities, both deterministic [8^ and stochastic [5]. 

As usual in problems of laser cooling, we will exten- 
sively use the fact that the time scales for the different 
processes are quite different. The first time scale is given 
by the time a photon needs to go through the sample. 
This time scale is so short that we can consider that the 
light follows immediately any change in the sample. A se- 
cond time scale is given by the atomic response time r~^ , 
which is in turn shorter than the third one related to the 
evolution of the external degrees of freedom (the atomic 
position and velocity). This inequality allows us to consi- 
der that the atomic internal state has always reached its 
steady state. 

In the following, we consider an atom in X with a ve- 
locity V, where we use capital letters to distinguish the 
position and velocity of this "real atom" from the coordi- 
nates {x,p) of a point in the phase space. The atoms are 
excited onaJ = 0— >J' = 1 transition by two counter- 
propagating laser beams with opposite circular polariza- 
tions. They also see a "bath" of photons which are spon- 
taneously emitted by all the other atoms in the sample. 
These photons propagate in all directions of the real 31? 
space and they have a broad spectrum. Although it is not 
always true, we consider here that the effect of the bath 
of spontaneous photons can be treated as a perturbation. 
Then, we split the force and the diffusion coefficient into 
two parts : one due to the trapping beams and the se- 
cond one due to multiple scattering : F = Ft + F„i and 



4.1 Effect of the trapping beams 



The trap consists in two counterpropagating beams 
with the same frequency ujl and opposite polarizations 
(Fig. [1^). The a~ beam has a wavevector k^ = k = ujl/c 
and has an intensity /_ {X, t) , while the a^ beam has 
a wavevector k^ = —k and an intensity I^{X,t). The 
frequency lol is slightly lower than the frequency wq of the 
atomic transition to insure an efficient Doppler cooling. 
The sample is submitted to an inhomogeneous magnetic 
field aligned on the propagation axis of the beams. For the 
sake of simplicity, we assume that the magnetic field varies 
linearly : B^ — bX. The Zeeman effect together with the 
cooling beams lead to a restoring force that gathers the 
atoms in the vicinity oi X = 0. 

In the atomic rest frame, the apparent frequencies of 
the trapping beams are Doppler shifted : a;± = wl ± kV. 
As the two beams have opposite circular polarizations, the 
relevant detunings are A± = A± S, where we have intro- 
duced the detuning A = uj^ — wq for an atom at rest with 
no magnetic field, and the sum 5 = kV -f ^sbX/h of the 
Doppler shift kV and the Zeeman shift fisBx/fi- As the 
two beams are circularly polarized, the excited state |0) 
plays no role in the present configuration. The coupling 
between the ground state \g) and the excited states |±) 
is described by the Rabi frequencies |/2±| = F^J I±/2Is, 
where Is is the saturation intensity. It is then possible to 
find exactly the steady state of the density matrix, by sol- 
ving the master equation. We are not interested here in 
the explicit form of the populations or of the coherences, 
which are given in a previous article [6]. The only im- 
portant point is that the stationary populations n± in 
the excited states |±) can be written as a function of the 
relevant parameters Z\, 5, /+ and /_. Note that these sta- 
tionary populations do not depend on the position or the 
velocity of the atom explicitly, but only through (5, I^ and 
/_. 

The key point is that the simplicity of the atomic 
structure means that 7± = Fn± are the emission rates 
of spontaneous a^ polarized photons. The exact expres- 
sions of n± can be found in [5]. Again because of the 
simple J = — ^ J' = 1 transition, a photon with the 
same polarization has to be absorbed before each spon- 
taneous emission process. As a consequence, we obtain a 
very simple expression for the radiative force due to the 
trapping beams, together with the expression for the dif- 
fusion coefficient due to the trapping beams : 



Ft = hk{-f^ -7+) 
A = ^/i'fc'(7-+7+) 



(8a) 
(8b) 



where a contribution 2fc^/5 in Dt comes from the mean va- 
lue of the square of the projected component of the spon- 
taneous wavevector. To evaluate this average, we consider 
the diagram of emission of a circular electric dipole (3/10 
in the plane and 2/5 in the orthogonal direction). 
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4.2 Effect of multiple scattering 

The contribution of multiple scattering is trickier to 
evaluate. We have to consider our atom in a bath of inco- 
herent photons which propagate in all directions. Because 
of their physical origin, photons do not differ in spectrum 
wherever the scatterer is and whatever their direction of 
propagation is. Therefore, the rate of absorption of pho- 
tons going in a given direction is simply proportional to 
the flux of photons that travel in this direction. The contri- 
bution to that flux of photons scattered just once is easy 
to evaluate. The flux of photons scattered more than once 
is more difficult to estimate, but it is in general much 
smaller than the previous one, and we neglect it in the 
following. The flux if'+(resp. 'F^) of photons, travelling to 
the right (resp. to the left) and impinging on an atom in 
X, is just half of the total number of photons scattered 
by the atoms on the left (resp. on the right) of X. 



1 



X 



^+{X,t)^-J dx' J dpp{x',p,t) (7- +7+) (9a) 
<P^{X,t) = -J dx' J dpp{x',p,t){j^+j+){9h) 

In these expressions, we do not differentiate the polariza- 
tion of the scattered photons. We just sum up the contri- 
bution 7_|_ of the mainly a^ polarized photons and the 
contribution 7_ of the mainly a^ polarized photons. It 
is interesting to note that the sum 'F^ + <F^ of these two 
fluxes does not depend on X. We shall see in the next 
section that these fluxes have much simpler expressions. 

To evaluate the force exerted by the scattered photons 
on the atoms, we have to know the fraction of photons 
which are absorbed and the momentum carried by each 
photon. This momentum is smaller than hk, because we 
have to keep here only the on-axis component for photons 
travelling in all directions of the real 3D space. Conside- 
ring the emission diagram of a circular dipole, we obtain 
(ka) = 9fc/16 and (kl) = 2k'^/5. This leads to the follo- 
wing expressions for Fm and Dm '■ 



9 

Id 
Dm^aR^r?k\F^+F+) 



(10a) 
(10b) 



where we have introduced re-absorption cross-section (Tr. 
In the limit of very low intensities (i? ^ 7^), we can 
give an estimate of ur. In this case, it is well-known that 
for a 2-level atom, the emission spectrum is mainly elastic 
P^ . This result holds for the J = ^ J' = 1 transition 
considered here and the scattered photons have the same 
frequency wl as the trapping lasers. To be rigorous, Dop- 
pler broadening should be taken into account, as the scat- 
ters are moving and the emitted photons do not propagate 
in the direction of the trapping beams. However, we ne- 
glect this broadening here, because it is smaller than both 
the detuning and the natural linewidth. We also neglect 
the Zeeman shift and the Doppler shift for the absorber. 
The absorption cross-section is thus reduced by a factor 



r'^ / (4Z\'^ + 1^^) with respect to the cross-section at reso- 
nance (To = 3A^/27r . In the very low intensity limit, we 
have : 

For higher but still modest intensities, a new contri- 
bution shows up : the blue sideband of the MoUow triplet 
|20) , which excites resonantly the atomic transition, grows 
up as ff I (4Z\^ + r"^) . As soon as J7 > T, this new reso- 
nant contribution dominates, even in the low saturation 
regime (1 > f^V (4/\2 + 7-2)). 

For even higher intensities, as we know the exact steady 
state of the density matrix, we can, in principle, compute 
the spectrum of the scattered light and the absorption 
cross-section, following what is usually done for 2-level 
atoms (llj . However, contrary to the previous publications 
[11113) , we do not evaluate or for a 2-level atom, to remain 
consistent with the model ofaJ = 0— )-J' = l transition 
where at least one of the excited state sub-levels does not 
interact with the trapping beams. As the exact calculation 
is quite heavy, we shall not make here the full calculation 
in the general case, but we shall restrict ourselves to the 
simpler situation considered in section [51 



5 Propagation of the trapping beams 

In this section, we derive the equation of propagation 
for the trapping beams assuming that we know the atomic 
density p {x, p, t). We are not interested in the phase of the 
laser beams, but just in their intensities. We have seen in 
the previous section that the rates 7-1- = rn± are the 
absorption rates of a polarized photons for an atom in 
X moving with velocity V. To know how many photons 
are absorbed between x and x + dx, we just have to sum 
the contributions of all atoms in X = x with all possible 
velocities V = p/m. These photons are taken from the 
trapping beam with the appropriate polarization, so we 
get : 



di+ 

dx 

di_ 

dx 



= +hojL I dpp{x,p,t)^+ (11a) 

-CXD 

= -hojL I dpp{x,p,t) 'f^ (lib) 



The plus sign in (|llap comes from the fact that the a^ po- 
larized beam propagates backwards. So /+ {x,t) increases 
with X, while /_ (x, t) decreases. 

Plugging (|lll) in the fluxes © and performing formally 
the integration on positions, we get : 



F+ {x, t) = 



/+ {X, t) - /_ (.T, t) - 1+ (-00, t) + /_ (-00, t) 

2hujL 

1+ (CX), t) - /_ (00, t) - 1+ {X, t) + /_ (.T, t) 

2hLJL 



The intensities /_ (—00, t) and /+ (00, t) are the incoming 
laser intensities, which are assumed to be constant. As 
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soon as we know the atomic phase space density p {x,p, t), 
we can solve the equations (fTTj) to get the spatial evolu- 
tion of the intensities. However, the evolution of the phase 
space density depends on the laser intensities through the 
radiative forces and the diffusion coefficients. Thus, the 
equations ([7]) and pT|) form a set of coupled equations 
which have to be solved together. 

It is important to note that Eqs pT|) introduce in the 
force a term which, in general, is not present in other 
systems described by the VFP equation. This difference 
has to be taken into account if we want to apply results 
obtained for plasmas to the dynamics of MOTs. Indeed, 
in plasmas, the particles do not act, by definition, on the 
thermal bath. On the contrary, in MOTs, atoms act on 
the beams (the equivalent of the thermal bath, as seen 
above), through the absorption and the so-called shadow 
effect. Thus it leads to another indirect dependence of the 
force on the density. 



6 Approximation for the atomic response 

In the general case, the equations ([7]) and (ITT|) descri- 
bing the MOT are highly non linear, and have no simple 
solution. Numerical simulations could give some solutions, 
but they could not give any physical insight of the phe- 
nomenon. However, starting from this general model, we 
can now consider some approximations allowing going fur- 
ther in the understanding of the MOT dynamics. Through 
these approximations, we aim to determine the expres- 
sions ([5]) of the trapping beam contribution and ([TUl) of 
the multiple scattering contribution to the mean force and 
diffusion coefficient of the MOT. 



6.1 Trapping forces 

Concerning the trapping beam components, Eqs ([8]) 
shows that we need essentially to explicit the scattering 
rates 7± = F n±, introduced as the emission rates of 
spontaneous a^ polarized photons, and also as the ab- 
sorption rates of the corresponding trapping photons. In 
the case ofaJ = 0— >-J' = l transition, these rates are 
exactly computable, as in [B]. However, the obtained ex- 
pressions are very heavy and difficult to interpret. To have 
a better physical insight, we will simplify these expressions 
through approximations valid in our domain of interest. 

The usual approach consists in considering the low sa- 
turation limit. In this case, each wave acts on the atoms 
independently, and the forces and the diffusion coefficient 
can be computed. But experiments are seldom realized 
in this low intensity limit. For slightly higher intensities, 
the calculation of the next order of perturbation (fourth 
order in field) must be considered. But this calculation 
is not straightforward. Indeed, the third order of pertur- 
bation builds up Zeeman coherence between the excited 
sublevels, which lead to non-trivial modifications of the 
populations at the fourth order. The calculation can be 



done, but the cross-terms do not allow any simple inter- 
pretation. Moreover, to evaluate the effect of multiple scat- 
tering, we need to know the absorption cross section for 
the bath of photons that have been scattered elsewhere. It 
requires one more step of perturbation by an extra weak 
field, which is not easy to implement. 

This first approximation is usually followed by another 
approximation, namely the low velocity limit, where the 
total trapping force is linearized in V. It naturally splits 
in a term due to the intensity imbalance and in a friction 
term. In the same way, the gradient of the magnetic field 
introduces a restoring force, linear in X. The condition 
of validity for this approximation is \S\ <c; max {\A\ , /^), 
where S is the sum of the Doppler and Zeeman shifts, as 
introduced in section 3] 

Let us consider here another approach, where the ap- 
proximation S is applied before going to the low saturation 
limit. The total shift 5 is considered as a perturbation, that 
is expanded to first order only. At a first glance, this choice 
seems quite surprising, as the "perturbation" is diagonal 
in the natural basis {\g) , |-|-) , |— )}. However, the good ba- 
sis to do the calculation is not the natural one, and then 
the shift S leads to off-diagonal terms. 

At zeroth order, the case i5 = corresponds to an atom 
interacting with a field which has both a^ components. 
It is thus natural to introduce the coupled \C) and non- 
coupled \N) states : 



\C) = 
\N) = 



[2: 



f2-\-) 



n*_ i-f) - nx I-) 



n 



where J7± are the generalized (C-number) Rabi frequen- 
cies that include the complex phase of the fields and fl = 

i7+ 1 -I- I Q- 1 describes the coupling between the ground 



state and the coupled one. Then the problem reduces to 
that of a 2-level atom interacting with a single laser field. 
The steady state population of the coupled state is sim- 
ply : 

(0) ^ n-" 

^ 4Z\2 + 2f22 + 7^2 
which gives the value of the scattering rates : 



(0) 

7± 



r ■ 



n\ 



4Zi2 



2r22 + r2 



The radiative force due to the beam imbalance and the 
coefficient of diffusion follows : 



^(0) ^ 



hkr 



4Z\2 + 2r?2 + r2 



(l^-l'-l^+P 



/?2 



* ~ 5 4Z\2 + 2122 + r2 



It is interesting to note that, in these expressions, we get 
the same result if we replace, in the expressions obtai- 
ned in the low saturation regime, the usual denominator 

(4Z\2 + ^2) by (4Z\2 + 2r22 -1- r2). 
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At first order, when 6 ^ 0, we introduce the same 
coupled and non-coupled states as above, and a first order 
perturbation in S gives the friction and restoring forces. 
The diagonal (in the natural basis) part of the hamiltonian 
Ha writes : 

HA = -iA+\+){+\+A^\-){-\) 
= -An\N){N\- Ac\C){C\ 

-24 i^+^- \N) {c\ + 12; r?* \C) {N\) 



with 



f2^ 

Ac = A + S 
An = A-S 



/22 



The slight change from A to Ac introduces a first order 
correction in the population of the excited coupled state. 



n. 



(1) 



(I 



AS \f2\^ -\f2: 



(4zi2 + 2122 ^ p2Y 



(13) 



The off-diagonal perturbation terms induce, at the first 
order, coherences between the non-coupled state and the 
two other states. But only the coherences in the excited 
state are needed to evaluate the populations n± 



r(l) 



r(l) 



= 



-2(5/2+l2_ 



4r 



4Z\2 + 2122 _^ p2 4/-^ _ j (2/^2 



N) = (N 



r(l) 



C 



(15) 



The correction in the scattering rates is deduced from the 
expressions of the population (fT5)) and of the coherences 
(HID and dH]) : 



7± 



T8A5r 



(4zi2 + 2/22 + p2y 



m\ 



2|/2_|^|i7+|^ 



4/^2 - 122 



i6r2z\2 + (2r2 + r22)'^ 



The evaluation of the friction force is then straightfor- 
ward : 



F, 



(1) 



hkrsAS 



(4Z\2 + 2/22 + 1^2)^ 



n^ +4\[2\^\[2:'^ 



4/^2 _ 1^2 



lenA^ + {2n + n'^y 



D 



(1) 



7 



ti^k^r- 



A5{\Q\' -\fi^ 



(4Z\2 + 2/22 ^ p2Y 



It is easy to check that the force calculated in this way 
is rather a friction force, as soon as Z\ < 0. The first 
term of the parenthesis, fP , together with the prefactor. 



corresponds to what is guessed from the unsaturated ex- 
pression : the original denominator (4Z\2 -\- f^^ becomes 
(4Z\2 -I- 2/22 _|_ p2'\^ However, a second term is needed to 
take into account properly the cross-saturation effects. 
This term is often forgotten by the authors, as for ins- 
tance in [S] , whereas it can modify significantly the spring 
constant of the trap, as it is of the same order of magnitude 
as /22. In particular, the minus sign means that, depen- 
ding on the parameters, the spring constant is increased 
or decreased significantly. On the contrary, the correction 
D\ to the coefficient of diffusion is usually very small. 



because |/2_ 



1/24 



6.2 Multiple scattering 



To evaluate the effect of multiple scattering, we ne- 
glect, as usual, the Doppler broadening due to the motion 
of the emitter and the Doppler and Zeeman shifts for the 
re-absorber. Thus, we consider that 5 vanishes. Eqs ([T0| 
shows that we have to evaluate the re-absorption cross- 
section gr, taking into account the fluorescence spectrum 
and the absorption spectrum, which are affected by the 
intense trapping field. 

When the laser field is intense enough, the dressed 
atom in the secular limit allows a simple evaluation of 
the fluorescence spectrum and of the re-absorption cross- 
section |11I19J . However, instead of considering a 2-level 
atom, we remain consistent with the model of a J = — > 
J' = 1 transition where we have three sublevels in the 
excited state. We have also to consider that the polari- 
zation of the photon scattered by an atom somewhere in 
the sample will not match exactly the polarization of the 
local field. For the sake of simplicity, let us consider here 
the case where the trapping beams have the same inten- 
sity, resulting in a linear polarization for the total trap- 
ping field. The consequence of the polarization mismatch 
is that, on average, one half of the scattered photons has 
the polarization of the local field while the other half has 
the orthogonal polarization. The detailed calculation is 
done in the appendix. The energy levels of the dressed 
atom are sketched in Fig. [3] In Fig. [3^, we have eviden- 
ced the fluorescence transitions that lead to the MoUow 
triplet. On the other hand, the absorption lines are drawn 
on Fig. [3]3, considering all the possible polarizations for 
the incident photon. Schematically, we have to consider 
the overlap of the four components of the Mollow triplet 
with the four absorption lines of the 3- level atom, as repre- 
sented in Fig. |3t. Some care has to be taken to compare 
the various contributions in the different regimes. Let us 
now examine the four main situations : 

(i) /2, |Zi| ^ r. In general, for large detunings and in- 
tensities and large light-shifts {5r = \ (%/Z\2 -I- /22 - |Zi|) > 
r) , the resonant term of the 2-level system dominates the 
sum, and all the other terms should be dropped to re- 
main consistent with the secular approximation done be- 
fore. The resulting re-absorption cross-section has thus the 
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Figure 3. Sketch of the energy levels of the dressed-atom. 
a) Spontaneous emission lines leading to the Mollow tri- 
plet, b) Absorption lines, c) Overlap of the fluorescence 
spectrum (solid line) with the absorption lines of the 2- 
level system (dotted lines) and with those due to the third 
level (dashed lines). 



following simple expression : 



Z\2|22 



<7R = 



0£ 

4 4Z\4 + 8Z\2/22 + 3124 



(16) 



This result is consistent, within a factor of 2, with the 
hmit values given for /2 > |Z\| > /" and |zi| > i7 > T in 
the appendix of Ref . ^Jj . 

(ii) \A\ :s> n. The contribution due to the non-coupled 
level is never exactly resonant. However, in this case, the 
light-shift, 6r — i7^/4 |Z\|, can be arbitrarily small, leading 
to an almost resonant behaviour ii Sr ^ F. In the case 
where |Z\| » i? > T and F > Q'^/4\A\, we take into 
account the quasi resonant contribution of the third level : 



CTR 






(17) 



Note that ^7^ is not the hmit of (HH) when \A\ > il, 
although it has the same form. 

(iii) Q <^ r. When the intensity is very small, the 
resonant or quasi-resonant contributions go to zero and 
the leading term corresponds to the reabsorption of the 
elastic component : 






(18) 



(iv) \A\ ^ r. For small detunings, the reabsorption 
by a 2-level atom goes to zero because the transition is 
strongly saturated. Then, the contribution of the third 
level, although it is non-resonant, has to be taken into 
account. We have : 



Ci? = Co 



25 /^ 
48 7?2 



(19) 



The above results allow us to write in the different re- 
gimes the contribution of the multiple scattering to the 
MOT dynamics. These expressions are rather simple as 
compared to those found in the literature. It is interesting 
to note that an. is always larger than the absorption cross- 
section for the trapping beams. This is very often assumed 
in previous papers, as stated in [13,. We have demonstra- 
ted here that, in the secular regime for three level atoms, 
it is always true. The above expressions of an, together 
with the expression of the trapping force derived in the 
previous paragraph, allow us to write a full model of the 
MOT, which can be used as a basis for future analyses. 



7 Conclusion 

In this paper, we built a model for a ID Magneto- 
Optical Trap as general as possible. We show that such a 
trap is described by a Vlasov-Fokker-Planck equation with 
a second relaxation term and a source term, both origina- 
ting in the bath of hot atoms of the atomic vapor. This 
VFP equation is coupled to a set of two differential equa- 
tions describing the beam propagation in the cold atoms. 
This system could be considered as relatively similar to 
plasmas, where the role of the thermal bath is played by 
the trapping beams. However, it appears that the MOT 
differs from plasmas on two important points : the second 
thermal bath, formed by the hot atoms, induces new inter- 
actions as compared to plasmas ; the trapping beams are 
not a "thermal" bath, as the atoms act on them through 
the absorption. 

In the last section of this paper, we derive in a more 
detailed way the equations established in the previous sec- 
tions, and compare the results with those found in the lit- 
terature. We show that important correction terms must 
be taken into account in the evaluation of the trapping 
forces, and we establish the expression of the multiple 
scattering in different situations. 

These equations aimed to be a basis for future de- 
velopments. They should contribute to obtain a better 
agreement between the experimental observations and the 
theoretical predictions. Let us remember that the previous 
existing models were only in qualitative agreement with 
the experimental observations, in particular for the situa- 
tions out of equilibrium. Using this model to study the 
instabilities in the cold atoms should give a better insight 
of the mechanism in which they originate. Thus a natural 
continuation of the present work would be to extend this 
model to the 3D traps. 



8 Appendix 

In this appendix, we present the detailed calculations 
of the re-absorption cross section, resulting from the over- 
lap of the fluorescence spectrum of an atom somewhere in 
the sample with the absorption spectrum of another atom 
elsewhere. In the J = — > J' = 1 transition conside- 
red here, we have three sublevels in the excited state, and 
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the polarization of the scattered photon will not match 
exactly the polarization of the local field. For the sake of 
simplicity, let us consider here the case where the trap- 
ping beams have the same intensity. In real experiments, 
it cannot be true everywhere in the sample because of 
the shadow effect, but the deviation remains on the order 
of 10%. In this simple case, the resulting polarization of 
the total trapping field is everywhere linear, but its direc- 
tion follows an helix. The laser wavelength is the scale for 
changing the relative phase between /?+ and f2_ , because 
the two beams are contrapropagating. As soon as we are 
not interested in what happens at length scales smaller 
than the laser wavelength, we can estimate that one half 
of the photons in the bath will interact with the transition 
between the ground state and the coupled one, while the 
other half of them can excite the atom to the non-coupled 
level. The total re- absorption cross-section, an, is thus the 
average of the usual cross-section, cr)j , calculated with a 
2-level atom [TT] and of a new contribution, a]j , coming 
from the presence of an excited level which is not cou- 
pled to the trapping field. The empty non-coupled state 
allows a strong absorption of the scattered light, while the 
absorption by the coupled state can be saturated by the 
trapping field. 

As usual, when one works with the dressed atom mo- 
del, the calculations are much simpler in the secular ap- 
proximation. This approximation is valid as soon as the 
energy splitting is larger than the typical width of the 
levels, which is expressed by the condition : 



fluorescence spectrum is given by 



S{u) = 



C4 + s4 



-5{uj -ujl) + 



2c^s 



,2^2^ 



c* + s^ 



-L{LLi-ujL,rp) 



+c''s^L{uj-ujL-nG,r,) 



where Hq = y/A^ + JP is the generalized Rabi frequency. 
The transformation from the bare basis to the dressed 
one is given by c = cos (9), s = sin (9), with the angle 
9 defined by tg {29) = — ^. 5 {ui ~ ujl) is the Dirac delta 
function, and L{uj,r) is the normalized Lorentzian : 



L{io,r) = 



r 



1 



TT LO^+n 



and the relaxation rates for the populations and the co- 
herences are : 



Tp = r (1 - 2c'^s^) 



r 

2" 



/; = - {i + 2c^s^) 



(22) 
(23) 



The first term in S {ut) corresponds to the elastic scatte- 
ring, while the three other terms are inelastic components. 
The last two terms are the sidebands of the MoUow tri- 
plet, centered in ujl ± ^Gi which are proportional to the 
intensity when \A\ ^ f2. 

The absorption spectra are given by : 



V^+^:>r 



(20) 



Please note that this relation does not require that both 
the detuning and the intensity are large, and some interes- 
ting limits can be examined with either i? < _r or |Z\| < /^, 
in the secular limit. First, we calculate the emission spec- 
tral density, S{lu), for an atom illuminated with the two 
trapping beams. Then, we evaluate the absorption spectra 
of another atom, also interacting with the trapping field. 
Two contributions appear in the absorption process : on 
one hand, the coupled excited level can still absorb light, 
leading to a contribution which is the one of a 2-level atom, 
<t\ (w), and on the other hand, the non-coupled level has 

to contribute with tj]^ (w). Finally, we evaluate the total 
re-absorption cross-section, or, with : 



1 



(2) 



.(3) 



^fl = 2 r« 



^f = I <^^A (w) ^ (w) '^^ 



(21) 



,(3) 



a\ (lo) S (w) du) 



All these calculations are done using the formalism of the 
dressed atom, in the secular limit (I^Ul) . The normalized 



4^^ (-) 



4^' (-) 



ttF c^ - s 



= ctq- 



2 c4 
ttF c^ - s' 



-c L{uj- UJL- ^G.Tc) 



-CTo- 



2 d^ + s' 



s'^L{uj - UJL + f^cTc) 



Oq—C L {lU - UlL - f^G +Sr,r2) 
+ CT0-^S L{UJ~LUL+ Sr, A) 



where Sr is the light-shift, and Fi^ are the relaxation rates 
for the coherences between the non-coupled level and the 
dressed states. 



Sr---(A + yf^^V^ 
^2 = f(l + ^^) 



(24) 
(25) 
(26) 



The absorption spectrum of the 2-level part consists in 
two lines centred in wl ± ^Iq , which means that they are 
resonantly excited by the two sidebands of the MoUow 
triplet. The two absorption lines due to the third level 
are centred in w/, -I- Qq ~ Sr and in uj^ — Sr, then the 
excitation by the fluorescence light is not resonant if Sr ^ 
r. Because of the secular approximation we made for the 
calculation, it is incorrect to write a^ as the sum of 16 
terms, coming from the overlap of the four components of 
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the fluorescence with the four absorption lines. We have 
to consider separately the different cases. 

In general, for large detunings and intensities (i7, |Z\| 3> 
r) and large light-shifts {Sr ^ -T), the resonant term of 
the 2-level system dominates the sum, and one should drop 
all the other terms to be consistent with the secular ap- 
proximation done before. The resulting cross-section has 
thus the following simple expression : 






A^n^ 



8Z\2i72 + 3l24 



(27) 



As already mentioned, the contribution due to the non- 
coupled level is never exactly resonant. However, if |Z\| 3> 
f2, the light-shift, Sr ~ il'^ /A \A\, can be arbitrarily small, 
leading to an almost resonant behaviour ii 6r '^ F . In the 
case where |Z\| > .!7 > T and F > fp- jA \A\, we take into 
account the quasi resonant contribution of the third level : 



^R 



16 A^ 



1 



16 A^ r2 



/?4 + 16 A^ n 



(28) 
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